A simple system with two temperatures 
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Abstract 

We study the stationary nonequilibrium regime which settles in when two single-spin 
paramagnets each in contact with its own thermal bath are coupled. The response vs. 
correlation plot exhibits some features of aging systems, in particular the existence, 
in some regimes, of effective temperatures. 
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A notion of effective temperature for non-equilibrium systems, first proposed 
by Hohenberg and Shraiman [1] in the context of turbulent flow, and later 
developed, among others, in ref. [2], is associated to the following Gedanken- 
experiment. A small system (called the thermometer) is coupled both to a 
heat reservoir at a temperature and to the system S whose effective tem- 
perature T eff must be determined. As a consequence, heat flows either from 
or to the reservoir at temperature 0. The effective temperature T e g is defined 
as the temperature at which this heat flow vanishes. In aging systems, as well 
as in glassy systems kept out of equilibrium by gentle stirring, it is expected 
that T e ff depends on the characteristic time scale of the thermomenter, and 
is related to the violation of the fluctuation-dissipation theorem. In order to 
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investigate this phenomenology, one has to develop a better grasp of the sit- 
uation in which a thermodynamical system is in contact with heat reservoirs 
at two different temperatures. 

Two-temperature systems naturally appear when considering systems with 
annealed degrees of freedom, like magnetic systems with evolving interactions 
[3] or diffusing particules on an evolving network [4], where spins or particules 
have a temperature different from the one of the bonds. We report here the 
results for what is possibly the simplest two-temperature system: two coupled 
single-spin paramagnets in contact with different thermal baths. When the 
temperatures of the baths are equal, the system reaches a trivial equilibrium 
state. When the temperatures are different, the system reaches a stationary 
nonequilibrium regime where energy flows through the system from the high- 
to the low-temperature thermal bath. We analyze this nonequilibrium state by 
means of the susceptibility vs. correlation plot [5]. At equilibrium, this plot is 
a straight line of slope —1/T, as a consequence of the fluctuation-dissipation 
theorem. In our case, the plot for the high-temperature paramagnet differs 
very slightly from the equilibrium one, but its effective temperature [2] is 
smaller than the one of its bath, in contrast with the more common situation. 
On the other hand, the low-temperature paramagnet can fall deeply out of 
equilibrium, and, depending on the parameters of the system, it can violate 
the fluctuation-dissipation theorem in a variety of ways, and are reminiscent, 
in some limit, of the corresponding plot in coarsening systems [6]. The ways in 
which the relation of the cross-response to the cross-correlation deviate from 
the equilibrium one are also interesting. 

We consider two classical single-spin paramagnets coupled via a bilinear spin- 
spin interaction. The Hamiltonian of the system is given by 

n = -\ E nSf + aSA (i) 

i=l,2 

where the "spin" variables Si can take any real value. The parameters r\ set 
the response time scale of the paramagnets and a represents the strength of the 
coupling. Stability requires a 2 < t\Ti- We shall consider a as a small parameter 
in the following. We recall that the equilibrium correlation Ci(t) = (Si(t)Si(0)) 
and response Ri(t) = 6 (Si(t)) /Shi(0) of an isolated paramagnet (where hi is 
a conjugate field to Si) are respectively given by 

Ci(t) = -e- ri|t| ; R i (t)=9(t)e- r ' t . (2) 



The dynamics of the coupled paramagnets is described by a set of linear 
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Langevin equations: 



dtSi = ~ + fh(t), 



(3) 



where 77, (i = 1,2) is a thermal noise, at temperature Tj, with zero mean 
and variance: (r)i(t)r)i(t')) = 2Tj<5(t — t'). Expliciting the previous equation we 
obtain 



8^ = -r 1 S 1 + aS 2 + r )l (t), 
d t S 2 = -r 2 S 2 + aS t + r) 2 (t). 



(4) 



iFrom this equation we derive a system of eight equations for the correlation 
Cij(t,tf) = {Si(t)Sj(t)) and response Rij(t,t') = 5 (Si(t)} / '6hj(t') functions 
respectively. 

For the response functions one has the autonomous equations 



(dt + rjRuit,?) = aR 21 (t,t') + 5(t-t'), 

(d t + r 2 ) R 22 (t, t') = aR 12 (t, t') + 5(t - t'), 

(d t + ri )R 12 (t,t') = aR 22 (t,t'), 

(d t + r 2 )R 2l (t,t') = aR n (t,t'). 



(5) 



By appying, e.g., (d t + r 2 ) to the third of these equations, and substituting 
the second one, we obtain 



{dt + r 2 ){d t + n) R 12 (t, t') = a 2 R 12 (t, t') + aS(t - f). 



(6) 



It is easy to see that R 2 \ satisfies the same equation. Since both i? 12 and R 21 
satisfy the same boundary conditions, namely they vanish for t < t', we deduce 
that 



R 12 (t,t') = R 21 (t,t'), Vt,t'. 



(7) 



For the correlation functions, the equations involve the response functions: 



(dt + rjCnit,?) = aC 21 (t,t') + 2T 1 R 11 (t f ,t), 
(d t + r 2 )C 22 (t,t') = aC 12 (t,t') + 2T 2 R 22 (t',t), 
(d t + ri )C 12 (t,t') = aC 22 (t,t') + 2T 1 R 21 (t',t), 
(d t + r 2 )C 21 (t,t') = aC u (t,t') + 2T 2 R 12 (t',t). 



(8) 
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After a short transient, the systems enters a stationary regime, where Cij(t, t') = 
Cij(t — t') and Rij(t,t') = Rij(t — t'). It is then possible to solve the system. 
We define the integrated response x(t) by 

t 

X (t)= J dt>R(t>). (9) 
o 



We then have 



xn(t) = e(t) 



(r 2 - a_)(e 



-a-t 



-I) 



a^(a + — cc_) 



+ 



(r 2 - a+)(e~ a + t -1) 
a + (a + — a-) 



Tl(r 2 _ Q 2 ) + a 2 T2 ^_ Q _ |t| Tl (r 2 2 -4) + a 2 T 2 |t| 



(a+ — a 2 ) 



— a 2 ) 



(10) 
(11) 



and the corresponding ones obtained by exchanging the labels 1 and 2. We 
have introduced the following notation for the inverse characteristic times: 



ri + r 2 

a ± = — — ^ 



n — ^2) 2 + 4a 2 



(12) 



Since the autocorrelation Cu(t) is a monotonically decreasing function of \t\, 
we can invert equation (11) for positive times, and express the integrated 
response in terms of the correlation. For times longer than a^ 1 , when the 
fastest decreasing exponential in eq. (11) can be neglected, we obtain 

y (t)- d( ®(l Q2 ^1 d{t) il ^ 1 (13) 

ah i ; Ti I a 2 + r2 ( ri+7 * 2 ^ / Ti \ a 2 + ( r2 ~ a -)( ri + 7 ' 2 ) ' 
\ t 2 /Ti— i / \ t 2 /Ti— 1 



We can directly read off this equation the slope of the response vs. correlation 
plot: 



X(C) d X _ 1 
Ti dC Ti 



\ 



a 



2 (r 2 - a_)(ri + r 2 ) 
a + 



(14) 



T2/T1 - 1 J 



This corresponds to a straight line, ie., to a well-defined effective temperature 
T eff = Ti/X(C) [2]. This temperature goes from Ti to infinity as T 2 /Ti grows 
from 1 to infinity, showing that the difference from equilibrium behavior can 
be very strong in the system coupled to the colder bath. 
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Fig. 1. Response function xu vs - correlation function C\\ for T± = 1 and different 
values of the higher temperature T 2 . We have set r\ = r 2 = 1 and a = 0.1. Both 
Xu and Cii are normalized by the instantaneous value Cn(0). The lines correspond 
(from above to below) to T 2 = 1, 30, 100, 300, 1000, respectively. 
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Fig. 2. Response function xu vs. correlation function Cu for T± = 1 and different 
values of the higher temperature T 2 . We have set r\ = 10, r 2 = 1, and a = 0.1. Both 
Xu and Cu are normalized by the instantaneous value Cn(0). The lines correspond 
(from above to below) to T 2 = 1, 100, 300, 1000, 2000, respectively. It is interesting 
to remark that in this situation the response dies off long before the correlation, 
like in coarsening. 

On the other hand, when T 2 /T\ goes from 1 to 0, the departure from equilib- 
rium is proportional to a 2 and positive, i.e., the plot lies slightly above than 
the FDT line. Thus the high-temperature paramagnet exhibits an effective 
temperature slightly smaller than the one of the bath. 

The expressions for the cross-correlation and response are: 



Xl2 {t) = a6{t) 
C 12 (t) = a 



-Ct-t 



-1 



a_(a + — a_) 
T l r 2 + T 2 r l f e~ a -^ 



a + (a + — a-) 

e -a+\t\ \ 



(15) 



a\ — a 2 _ 
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(16) 

and the corresponding ones obtained by exchanging the labels 1 and 2. For 
times longer than a^ 1 , when the fastest decreasing exponential in eq. (16) can 
be neglected, the cross-susceptibility has a simple expression in terms of the 
cross-correlation: 



_ sign t(T!-T 2 ) , Q _ |t| _ e _ Q+|t| 
a\ — o?_ ^ 




n+r 2 ^ , a_(7i-T 2 ) \ 



riT 2 + r 2 Ti V riT 2 + r 2 Ti - a_(7i - T 2 ) 



Notice that the curves of the cross-response vs. the cross-correlations start 
from the same point (since Ci 2 (0) = C 2 i(0), while Xi2(0) = X2i(0) = 0) and 
end at the same point, because the symmetry of Rij(t) implies the equality 
of Xi2(t) and X2i(t)- Nevertheless the two effective temperatures are different, 
one being above and the other below the "average" temperature T = (r{T 2 + 
r 2 T 1 )/(r 1 + r 2 ). 
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Fig. 3. Cross-response function xn vs. correlation functions C7 12 and C 2 i for T\ = \ 
and different values of the higher temperature T 2 . We have set n = 10, r 2 = 1, and 
a = 0.1. Ci 2 , C 2 \ are normalized by the instantaneous value C'i 2 (0) = 6*21(0), and 
X12 by its asymptotic value for t — > 00. The lines correspond (from outside to inside) 
to T 2 = 10,5,2,1, respectively. The lines above correspond to C\ 2 , those below to 
C 2 \. Notice that the behavior of C\ 2 is not monotonic. 

Summarizing, we have shown that the behavior of a very simple system with 
two temperatures exhibits some of the characteristic features of aging systems, 
kept out of equilibrium by a stirring force. It is possible to extend the approach 
to study the details of the measurement of temperature in an aging system [7]. 

We thank Prof. Serge Galam for his interest in this work. L. P. acknowledges 
the support of a Chaire Joliot of the E.S. P.C.I. 



6 



References 



[1] P. C. Hohenberg and B.I. Shraiman, Physica 37, 109 (1989) . 

[2] L. F. Cugliandolo, J. Kurchan and L. Peliti, Phys. Rev E55, 3898 (1997). 

[3] H. Horner, Z. Phys. B57, 29, (1984); ibid. 39, (1984). A. E. Allahverdyan and 
D. B. Saakian, Phys. Rev. E58, R5201 (1998). 

[4] A. Fierro and A. Coniglio, in preparation. 

[5] L. F. Cugliandolo and J. Kurchan, J. Phys. A27, 5749 (1994). 

[6] L. F. Cugliandolo and D. S. Dean, J. Phys. A: Math. Gen. 28, 4215 (1995). 
C. Godreche and J.-M. Luck, J. Phys. A: Math. Gen. 30, 6245 (1997). A. Barrat, 
Phys. Rev. E57, 3629 (1998). 

[7] R. Exartier and L. Peliti, |cond-mat/9910412 



7 



